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Preface 


This booklet contains a collection of contributions to the meeting of the JEts and Disks at INAF 
(JEDI) group, which took place at the Capodimonte Observatory during 9-10 April 2015. Scope 
of the meeting was to bring together the JEDI researchers of the Italian Istituto Nazionale di 
Astrofisica (INAF) working in the field of circumstellar disks and jets in young stars, to discuss 
together the different agents affecting the structure and the evolution of disks, namely accre¬ 
tion, jets and winds. More information on the JEDI group and its activities can be found at 
http://www.oa-roma.inaf.it/irgroup/JEDI. The JEDI group is supported through funds from 
the PRIN-2012 “Disks, jets and the dawn of planets 
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Part 1: 


Early Protostellar Phases 



THE CH3CHO/HDCO RATIO AS A TOOL TO STUDY THE COMS FORMATION 


E. Bianchi^, C. Codella^, F. Fontani^, L. Testi^’^, R. Bachiller®, C. Ceccarelli^’®, B. Lefloch'^’®; ^INAF - Osservatorio Astrofisico 
diArcetri, Firenze, Italy, ^European Southern Observatory, Garching, Germany, ^IGN Observatorio Astronomico Nacional (IGN), 
Madrid, Spain, '^Univ. Grenoble Alpes, IPAG, Grenoble, France, ^CNRS, IPAG, Grenoble, France. 


Understanding the emergence of molecular complexity, 
from simple atoms and molecules to the large number of chemi¬ 
cal species observed in the Solar System, is a crucial astrophy s- 
ical topic. Complex Organic Molecules (COMs, > 6 atoms) 
have been detected in all the components of the star formation 
process, demonstrating the existence of efficient pathways to 
chemical complexity. During the low-mass star forming pro¬ 
cess, the primordial cloud evolves into a protostellar envelope, 
protoplanetary disk and planetary system. At the same time 
the interstellar medium evolves toward denser condensations 
and molecules can therefore be formed, destroyed or incorpo¬ 
rated at the varius stages so that the chemical composition of 
the gas becomes increasingly more complex. 

One of the simplest COMs is acetaldehyde (CH3CHO). 
It has been recently detected in the chemically rich outflow 
driven by the LI 157 mm protostar (Class 0), and associated 
with molecular clumpy cavities created by episodic events 
in a precessing jet. In particular, acetaldehyde has been re¬ 
vealed toward its brigthest bow-shock called B1. Interestingly 
CH3CHO shows the same spatial distribution of deuterated 
formaldehyde (HDCO), both following the young outflow cav¬ 
ity produced by the impact of the jet with the ambient medium 
[4,3]. These observations indicated that the two molecules are 
closely associated with the region enriched by species evapo¬ 
rated from icy grain mantles and released into the gas phase. 

In order to verify this result, we aim to extend the sample 
analyzing other three young stellar objects: SVS13 A & B, 
Cep E and HH212. We present here the preliminary results 
of the deuterated formaldehyde (HDCO) observations toward 
SVS13. 

This object is located in the NGC1333 star forming region 
at a distance of 235 pc [5]. It is associated with several young 
stellar objects and in particular with two main bright mm- 
sources called A and B. SVS13Ais toughtto be a Class I object 
while SVS13B should be younger (Class 0) as suggested by 
the occurrence of a collimated jet [1]. 

Spectra of SVS 13A were observed with the IRAM 30-m 
antenna, in the framework of the ASAI large program (PI: B. 
Lefloch & R. Bachiller), an unbiased spectral survey aimed to 
observe several targets sampling different stages of solar-type 
star formation. As a first step we searched for the HDCO 
rotational transitions in the spectra at 1mm, 2mm, 3mm, us¬ 
ing the GILD AS package'and the JPL^ database to identify 


the lines. We clearly detected four lines in the 1mm spectra 
one at 2mm, while the 3mm band did not present any HDCO 
transition. This could be due to a beam filling factor effect 
(HPBW 10”, 20”, 30” at 1, 2, 3 mm, respectively). If the 
emitting size is significantly smaller than the beam size then 
the emission could be diluted. Note that the HPBW at 1mm 
is smaller than the angular separation (<^14.5 arcsec) between 
the two main sources and samples only SVS 13A. The prelim¬ 
inary analysis of the 1mm spectra shows four spectral lines 
(up to Eu ”^40 K) with a peak velocity of 8 km/s and 
typical FWHM of 1-A km/s. Moreover some lines present an 
indication of a weak broader component that could be associ¬ 
ated to an outflow. In Fig 1 we show the lines profile and the 
gaussian fit. In order to determine the excitation conditions 
i.e. excitation temperature and column densities, we build 
rotational diagrams of the different line components. This 
approach assumes that the lines are optically thin, and Local 
Thermodynamic Equilibrium (LTE). In Fig 2 we show the ro¬ 
tational diagrams. The green points indicate the bulk of the 
central line narrow emission around the peak velocity while 
the blue and the red points correspond respectively to the blue 
and the red broader wings. We found rotational temperatures 
of 20 K and column densities of A 4 x 
for the central emission, A 3 x for the blue 

wings and A 5 x for the red one. Further 

analysis will be carried out for the detection of the acetalde¬ 
hyde (CH3CHO) in the same source and subsequently of both 
HDCO and CH3CHO in the other objects. 
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E. Bianchi - The CH 3 CHO/HDCO ratio as a tool to study the COMs formation. 




Figure 1: Deuterated formaldehyde (HDCO) detected spectral lines toward SVS13A. The lines were observed at 1 mm with 
the IRAM 30-m antenna, in the ASAI framework (PI: B. Lefloch & R. Bachiller). The spectra are reported in the main beam 
temperature scale. The peak velocity is ~+8 km/s, while the typical FWHM are 2-4 km/s. In red we show the gaussian fit of the 
spectral lines. 
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Figure 2: Rotational diagram for HDCO as observed toward SVS13A. The green points indicate the velocity range of <^5 km/s 
around the peak velocity. The blue (red) points correspond to the blue (red) wings. The black lines are for the best fits. The values 
of both rotational temperatures and coloumn densities are derived under the assumption of optically thin emission and LTE. For 
the red wings we assume the same rotational temperatures of the blue one. 
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HH212’S ANATOMY: STORY OF A PROTOSTELLAR SERIES 


C. Codella; INAF - Osservatorio Aslrofisico di Arcetri, Firenze, Italy. 


Magneto-hydrodynamical (MHD) driven jets play a cru¬ 
cial role during the earliest Class 0 phase, that is in the star 
and disk formation process itself. High-resolution MHD sim¬ 
ulations of protostellar collapse show that MHD outflows are 
an unavoidable outcome, and are able to eject 20%-50% of 
the infalling core gas before it reaches the central source [ 2 ]. 
The magnetic braking by the outflows and twisted B-fields is 
so efficient when the field and the spin axis are aligned that 
keplerian disks may be initially suppressed beyond 10 AU 
(e.g. [9]). On the other hand, large keplerian disks of 100- 
150 AU have recently been reported in three Class 0 sources 
[11,10,8]. Such large Keplerian disks might result from mis¬ 
aligned B ~ Q. configurations or a strong turbulence [5], but 
observational evidence is still lacking. A detailed observation 
and characterisation of a larger sample of Class 0 systems is 
thus essential to elucidate the disk formation process. 

HH212 is a bright symmetric bipolar jet from a Class 0 
source in Orion (e.g. [13]). Previous Sub-Millimeter Array 
(SMA) and Plateau de Bure Interferometer (PdBI) maps of the 
inner (< 900 AU) region revealed a bipolar SiO/CO microjet 
invisible in H 2 due to high extinction (e.g. [3,6,1]). The 
dust continuum recently imaged by ALMA suggests a disk of 
(maximum) radius ~ 0.3 AU [7]. All these findings make 
HH212 an ideal laboratory for investigating the interplay of 
infall, outflow and rotation in the protostellar stages. 

We summarise here the first results obtained using HH212 
observations performed during the ALMA Cycle 0 phase. The 
CO(3-2), C^'^0(3-2), SiO(8-7), C®'‘S(7-6), andCH30H(7i.7- 
61 , 6 ) lines at 345795.99 MHz, 337061.13 MHz, 347330.63 
MHz, 337396.69 MHz, and 335582.00 MHz respectively, were 
observed using spectral units of 0.43 km s“^ resolution. Fur¬ 
ther details can be found in [4]. 

Figure 1 compares the emission maps in the 850 /rm 
continuum, SiO(8-7), CO(3-2), C^'^0(3-2), C®'‘S(7-6), and 
CH 30 H( 7 i, 7 - 61 , 6 )- The combination of these tracers allows 
us to simultaneously image in a single ALMA-Band 7 spectral 
set-up different ingredients of the star formation process: (i) 
a pair of narrow SiO jets launched from the protostar (dust 
and CH 3 OH peaks), (ii) an extended flattened C^^O envelope 
around the outflow waist, and (iii) biconical CO and C®^S 
emission lobes surrounding the jet. 


Figure 2 shows that, in addition to infall motions (traced 
by the extended diamond shape of the position-velocity plot), 
a rotation signature is clearly seen in the form of two emission 
peaks at low velocities (LV; < 1.5 km s“^): one blueshifted 
to the west, and one redshifted to the east. This rotating 
LV C^^O emission is tracing the sides of the southern cavity 
carved by the outflow into the envelope, which rotates in the 
same sense [12]. At high velocities (HV; ~ 1.9-3.5 km s“^ 
from systemic): (i) the E-W velocity gradient is still present, 
and (ii) the emission is definitely more compact (<0.5 arcsec). 
We obtained the emission centroid positions in each channel 
from fits in the uv domain. While the centroids in the LV 
range fall on the southern cavity, the centroids in the HV 
range move to the equatorial plane, which indicates an inner 
rotating equatorial disk of radius ~ 0.2 arcsec (90 AU) around 
a 0.3±0.1 Mq protostar. 
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C. Codella - HH212’s Anatomy 
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Figure 1: The HH212 protostellar system as observed with ALMA-Band 7 [4]. Left panel: Contour plot of the continuum emission 
(turquoise contours) at 0.9 mm overlaid on top of the CO(3-2) emission (colour scale and black contours), integrated in the Vsys 
±17 km s“^ range. The ellipse shows the ALMA synthesised beam (HPBW). Labels indicate the MMl protostar. Middle panel: 
Contour plot of the C^^O(3-2) emission (turquoise) on top of the C®^S(7-6) emission (colour and black), both integrated in the 
14aB±6 km s“^ range. Right panel: Contour plot of the CFl 30 H( 7 i, 7 - 61 , 6 ) emission (turquoise, integrated in the ’I4j,s±6 km s“^ 
range) on top of the SiO(8-7) emission (colour and black, integrated between -23 and +15 km s“^ with respect to Vsya)- 
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Figure 2: Upper-left: Blue- and redshifted C^^O(3-2) emission in the low-velocity (LV) range (±0.6-1.5 km s“^ from Vsya), 
tracing the rotating outflow cavity. The tilted black cross indicates the SiO jet direction (PA = 22°) and the equatorial plane. 
Upper-right: same as the Upper-left panel for the high-velocity (FIV) velocity range (±1.9-3.5 km s“^ from systemic), tracing 
the rotating inner disk. Bottom-left: position-velocity cut of C^^O(3-2) perpendicular to the jet. Dashed lines mark Vaya and the 
continuum peak MML Bottom-right: distribution of the C^^O(3-2) centroid positions (from fits in the uv domain with la error 
bars) in the various velocity channels. Magenta is used for the systemic velocity. Red and blue datapoints denote the channels 
whose fitted centroids lie, within the error bars, on the equator. 
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DEUTERATED MOLECULES: A CHEMICAL EILTER EOR SHOCKS 


F. Fontani; INAF - Osservatorio Astrofisico di Arcetri, Firenze, Italy. 


The process of deuterium enrichment in molecules from 
the main reservoir HD is initiated in molecular clouds by 
exothermic gas-phase reactions, and it is thus boosted in cold 
and dense gas (e.g. [1]). For this reason, deuterated molecules 
have long been used as tracers of the coldest and densest por¬ 
tions of interstellar clouds. Although this general framework 
remains true, increasing observational evidence is showing that 
the story is not so simple: deuterium fractionation is boosted 
in cold and dense gas, but it does not uniquely trace cold 
and dense gas. In fact, a second very efficient way to form 
molecules rich in hydrogen and deuterium is on the surfaces of 
dust grains during the early cold pre-stellar phase, which are 
then released into the gas by thermal or non-thermal desorp¬ 
tion mechanisms. Therefore, their abundance can reach high 
values even in hot gas (see [2] for a review). High deuter¬ 
ated fractions (i.e. the abundance ratio of the isotopologue 
containin D and its counterpart containing H) of species that 
are efficiently formed on dust grains, like formaldehyde, am¬ 
monia, and methanol, are indeed measured in hot-cores and 
hot-corinos (e.g. [2]), where evaporation of grain mantles due 
to heating of a protostellar object is maximum. However, up 
to now, the deuterated fraction in shocks was not investigated. 

The first source in which this investigation has been per¬ 
formed is L1157-B1, the brightest bow-shock associated with 
the protostellar outflow LI 157 (the prototype of chemically 
active outflows) driven by a low-mass class-0 protostar [3]. To¬ 
wards LI 157-B1, several emission lines of deuterated molecules 
have been detected by our team for the first time (including 
HDCO and CH 2 DOH) with the IRAM-30m Telescope [4]. 
Based on chemical models [5], we suggested that these species 
were formed on grain mantles, and then released into the gas 
phase by the passage of the shock. The young age of the shock 

2000 yrs, [6]) is consistent with this picture, because in 
such a short lifetime the chemical reactions in the hot gas un¬ 
likely can modify significantly the chemical composition of 
the evaporated material (e.g. [7]). However, the angular reso¬ 


lution of the data published by [4] was not sufficient to firmly 
claim that the detected emission indeed arises from the shock 
impacting region. This scenario has been fully confirmed 
by observations obtained at high angular resolution with the 
Plateau de Bure Interferometer [8]: as expected the emission 
of HDCO perfectly delineates the region of interface between 
the fast jet and the slower ambient material (Fig. 1). CH 2 DOH 
emission is fainter and thus its emitting region is not as clearly 
delineated as for HDCO, but its integrated map is consistent 
with gas mainly ejected where the jet impinges the ambient 
material. The deuterated fraction HDCO/H 2 CO is 0.1 in 
the HDCO emitting region, an order of magnitude larger than 
the upper limit found in the surrounding material, probably 
dominated by warm-gas chemistry and less affected by grain 
evaporation. Our study represents the first clear evidence ever 
found of HDCO as shock tracer, proposing a new use of this 
molecule as "chemical filter" to trace recently evaporated ma¬ 
terial (more selective than its hydrogenated counterpart, as it 
can be seen from Fig. 1). Moreover, the study yields an in¬ 
direct but "clean" measurement of the deuteration of the ices 
covering the dust grains during the cold pre-protostellar phase. 
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F. Fontani - Deuterated molecules 
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Figure 1: (a): map of the intensity of the HDCO( 2 i 4 — lo,i) line (thick black contours) averaged between -0.7 and 3.4 km s“^, 
observed with the PdBI. Contours are 3, 5 and 7 times the Icr rms level in the map, i.e. 3.8 mjy beam“^. The grey scale 
represents the averaged emission of H 2 CO( 2 o ,2 — lo,i) observed with the PdBI by Benedettini et al. [9]: the levels range from 15 
(3 c7 rms level) to 190 mJy beam“^, in steps of 30 mJy beam“^. The big circle represents the PdBI primary beam 37") at the 
frequency of the HDCO line, and the ellipses in the bottom- and top-left corners show the PdBI synthesised beam for the HDCO 
and H 2 CO images, respectively, {by. SiO(2-l) (heat colour scale) and CH 3 CN( 8 - 7 ) (K=0-2, cyan contours) averaged emission, 
with the location of the CH 3 CN condensations identified by Codella et al. ([10], same contour levels as in their Fig. 1 are shown). 
To the left of the maps, we show the spectra of HDCO(2i,i — lo,i) and H 2 CO( 2 o ,2 — lo,i) integrated over the red, blue and 
green contours (from top to bottom) in panel (a). These delineate: the eastern wall ("E-wall”) of the cavity detected in HDCO, 
the arch-like structure ("arch”) detected in HDCO behind the head of the bow-shock, and the ’’head” of the bow-shock detected 
in H 2 CO but not in HDCO. The purple line in the HDCO spectra shown in the ’’E-wall” and ’’arch” panels represents the best 
Gaussian fit to the lines. From [ 8 ]. 
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SEARCH FOR ACCRETING PLANETS IN TRANSITIONAL DISCS 


J.M. Alcala^, E. Covino^, S. Antoniucci^, B. Nisini^, T. Gianinni^, C. Manara^; ^INAF - Osservatorio Astronomico di Capodi- 
monte, Napoli, Italy, ^INAF - Osservatorio Astronomico di Roma, Monte Porzio Catone, Italy, ^ESA/ESTEC, Noordwijk, The 
Netherlands. 


A key issue in the study of planet formation is to explain 
how optically thick accretion discs surrounding the youngest 
solar-mass stars, i.e. the classical T Tauri stars, evolve into 
optically thin debris discs around young main sequence stars, 
passing through the weak-T Tauri star phase [22]. Gener¬ 
ally, these discs persist for several million years [ 11 ], during 
which part of the material is accreted onto the star, part is lost 
via outflows and photoevaporation, and part condenses into 
centimeter-sized and larger bodies or planetesimals. An in- 
termdediate stage of T Tauri disc evolution is observationally 
identified with the so-called "transitional discs" (TDs), which 
are characterised by inner holes in their dust distribution. The 
analysis of the spectral energy distribution (SED) and interfer¬ 
ometric measurements at sub-millimeter wavelengths indicate 
inner cavity radii, Ri„, spanning/rom a few AUs to several 
tens of AU [9,2,16,10]. 

Important clues to the process of planet formation are ex¬ 
pected from observation of TDs, as they represent the best 
candidate places where forming planets can be found, located 
inside the circumstellar disc cavity and close to the inner rim 
of the disc [7]. Therefore, the estimates of Ri„ provide upper 
limits for the expected orbital radius of the planet, although de¬ 
pending on the SED modelling approach used [12,20]. There¬ 
fore, direct detection of the inner rim of TDs is also crucial to 
costrain the models. 

I. Accreting planets in TDs 

Recent measurements of high mass-accretion rates {Macc^ 
lO“®-lO“®M 0 yr“^; [9,15,1]) in several TDs (see Fig. 1) 
suggest that there may still be an important amount of residual 
circumstellar gas in the inner cavity. Interestingly, no corre¬ 
lation is found between Macc and the size of the inner cavity 
(Fig. 1, bottom panel), as would be expected for magneto- 
spheric accretion onto the star. This suggests that the density of 
the gas in the inner regions is independent of the mechanisms 
producing the cavity [15]. Therefore, the accretion activity 
has also been ascribed to mass accretion onto forming planets 
[14,18,24]. Indeed, since circumstellar discs have lifetimes of 
some million years, forming a 1-10 Mj„p planet during this 
time lapse requires that a circumplanetary disc has an average 
accretion rate of Macc> lO“®-lO“®M 0 yr“^. However, ex¬ 
cept for the candidates around LkCa 15 [13], HD100546 [17] 
and HD169142 [19,5], objects showing evidences of on-going 
planet formation have not been clearly detected yet. 

2. Detection of circumplanetary discs 

Recent calculations [24,8] of the SEDs of circumplanetary 
accretion discs have shown that a circumplanetary disc with 
Macc~ lO“®M 0 yr“^, around a 1 Mjup planet has a max¬ 
imum temperature of <^2000 K, and in the near-infrared it 


should be as bright as a late-type brown dwarf or as a 10 M 
planet . Hence, such a disc could be detected by using the 
high-contrast, high angular resolution and high sensitivity in¬ 
strument SPHERE @VLT. However, in order to identify un- 
ambiguosly circumplanetary accretion and distinguish it from 
brown dwarf companions or hot high-mass planets, we need 
to reveal accretion signatures such as Ha emission, and/or 
find indication from mid-infrared photometry of a flatter SED, 
caused by a warm circumplanetary disc, as compared to the 
SED of brown dwarfs or non-accreting planets [24]. Recent 
studies have stressed on the use of Ha emission as a powerful 
diagnostic to find both accreting stellar companions [ 6 ] and 
proto-Jupiters [23]. 

3. The project 

Detecting and characterising planetary-mass companions of 
young stellar objects (YSOs) with TDs and studying the dust 
distribution in the inner disc rims are primary observational 
steps for unveiling planet formation processes and for cos¬ 
training the models of late disc evolution. Therefore, within 
the framework of lEDI we have started a project to investigate 
YSOs with TDs using SPHERE® VLT and LMIRCam@LBTI, 
in close collaboration with the INAF-Padova group incharge of 
high-contrast, high spatial resolution imaging techniques. The 
main goal of this project is to look for planetary-mass com¬ 
panions around YSOs with TDs. The YSOs span a narrow 
range of ages, between 1 Myr and 3 Myr, and show evidence 
of mass-accretion, which makes them the best candidates to 
host forming planets. Planetary companions are expected to 
be located inside the cavities of the TDs and close to the inner 
edge of the disc cavity [7]. 
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Figure 1: Top: typical SED of a TD for the YSO Sz91 in 
Lupus III (black symbols). A model atmosphere of the same 
effective temperature as the star is overplotted (blue line). The 
’’valley” between 10 pm and 100 pm due to the lack of dust 
in the inner disc regions is clearly seen. The estimated size 
of the inner cavity is Ri„ = 65 ALf ([21]). Middle: UVB X- 
Shooter spectrum of Sz91 (red) showing the noticeable accre¬ 
tion signatures as UVB excess emission and Balmer jump with 
respect to a template discless star of the same spectral-type 
(green). A mass accretion rate Macc=L5x 10“® M 0 yr“^ 
was calculated by fitting (blue line) the template plus the emis¬ 
sion of a slab of hydrogen (black) (adapted from [1]). Bottom: 
Macc versus radius of the inner cavity of a sample of TDs 
(adapted from [15]). 
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HI DECREMENTS AND LINE PROFILES IN YOUNG STELLAR OBJECTS 


S. Antoniucci; INAF - Osservatorio Astronomico di Roma, Monte Porzio Catone, Italy. 


Hydrogen recombination lines (HI) are one of the most com¬ 
mon features observed in optical and near-infrared spectra of 
young stellar objects (YSOs) and represent a powerful tracer 
for the high-density gas located in the innermost circumstellar 
region of young stars. These lines are commonly used as a 
tracer of the mass accretion process [1,2], but they can also 
be effectively employed to probe the physical properties of 
the gas in the circumstellar structures [6,10,11]. Instruments 
like VLT/X-Shooter, with its broad wavelength coverage and 
moderate spectral resolution, allow us now to simultaneously 
study HI lines from the Balmer, Paschen, and Brackett se¬ 
ries, and to derive both the stellar and accretion properties of 
the young sources from the analysis of spectral features and 
accretion-related excess emission [1,9]. 

In our work [3] we have considered the X-Shooter spectra 
presented by [1], which refer to a sample of 36 low-mass 
Class II sources in the Lupus low-mass star-forming region [5] 
for which all the stellar and accretion properties are already 
available. Thanks to the good quality of the UVB and VIS 
spectra of this X-Shooter dataset, we have focused on the 
lines of the Balmer series, and in particular on the analysis 
of the series decrement. On the one hand, we have looked 
for correlations between the decrement shape, the observed 
line profiles, and the (stellar and accretion) properties of the 
objects; on the other hand we have compared the observed 
decrements to the predictions of standard emission models 
(case B [7], and the local line excitation computations by [8]). 
Our aim is to investigate what are the gas properties in the 
line emission region and whether these vary throughout the 
sample. 

I. Decrement shapes 

The extinction-corrected Balmer decrements of the targets rel¬ 
ative to the H/3 (H4) line and up to the H15 line are shown in 
Fig. 1. Although the decrement shapes vary with continuity 
within the sample, we have identified four major decrement 
types, which are reported in the panels of Fig. 1, using a rep¬ 
resentative source for each type. Only three objects present a 
type 1 decrement, which presents a curved shape outlining the 
lower part of the decrement ensemble. Type 2 decrements are 
the most common (20/36) and look as roughly straight lines in 
the plot. Type 3 decrements (7/36) have a shape characterised 
by a fairly wide bump and unlike type Is occupy the upper side 
of the decrement ensemble. Finally, in type 4 (6/36) the slope 
is steep for the first lines of the series and becomes flatter for 

N„p > 6. 

2. Decrements and source properties 

We note that the three sources displaying type 1 decrements 
have been previously identified as sub-luminous objects [1], 
most likely due to geometrical obscuration of the central region 
from an edge-on view of the disc. In such a geometry it is likely 


to have a residual amount of (differential) extinction that is not 
properly accounted for, which would effectively result in a 
type 1 shape decrement. 

We could not determine any clear correlation between the 
other decrement shapes and the source properties, except for 
a tentative connection with the mass accretion rate. Indeed, 
the strongest accretors in the sample seem to show similar 
decrement shapes: type 4 (mostly) and type 2. 

3. Decrements and line profiles 

Twelve targets present highly symmetric lines with a FWHM 
around 100 km/s {narrow symmetric line objects). Much larger 
profiles (FWHM up tp 350-400 km/s) that are still symmet¬ 
ric are observed in eight sources, while the remaining sources 
show more complex (double- or triple-peaked) profiles, as a re¬ 
sult of strong absorptions. The FWZI measured in the sample 
are in the range 400-1600 km s“^. The comparison between 
FWHM and FWZI indicate that most objects have lines that 
significantly differ from a Gaussian, having very extended 
wings, which may signal gas moving at very large velocities 
with respect to the star. The most evident connection between 
decrement shape and line profiles is that all sources with nar¬ 
row symmetric lines display a type 2 decrement, although not 
all type 2 decrements are associated with narrow symmetric 
line objects. Targets with type 4 decrements have in gen¬ 
eral very wide profiles, which suggests emission from a fairly 
extended circumstellar region in which large amounts of gas 
move at different velocities, a picture that is consistent with 
the connection between type 4 decrements and higher mass 
accretion rates. 

4. Decrements and models 

Case B emission is able to reproduce well type 2 and type 
3 decrements assuming large densities (> 10^° cm“®) and 
low temperatures (1000-2000 K) [4], but in order to satisfy 
the assumptions of case B (optically thin lines) under these 
conditions a very strong source of ionising photons is required, 
which is hardly found in circumstellar regions of typical T 
Tauri stars. 

Models by [7] can manage optically thick lines, although 
in a local environment, i.e. assuming a gas with a single tem¬ 
perature and density, a velocity gradient, and an ionisation 
rate. These models reproduce well sources with type 4 decre¬ 
ments, indicating (total) densities around 10^°-10^^ cm~®) 
and temperatures around 10000 K. This suggests that in this 
sources the bulk of the emitting gas is well described in terms 
of a single temperature and density. The most common type 
2 decrements are reproduced less well at somewhat lower gas 
densities (10®-10^° cm~^), while the temperature remains dif¬ 
ficult to constrain. Type 3 decrements remain the most difficult 
to reproduce in the framework of these models. 

A more complex modelling of the emitting region (includ¬ 
ing a temperature and density profile) is probably required to 
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Figure 1: Balmer decrements of four sources of the Lupus sample. Each target represents one of the four main decrement shapes 
we have identified, which are plotted in red against the ensemble of all decrements of the sample (grey curves). Adapted from [2]. 


fully reproduce the observed decrements and derive the prop¬ 
erties of the emitting gas. 
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A STUDY OF ACCRETION IN THE L1615/L1616 COMETARY CLOUD 


K. Biazzo^, J. M. Alcala^, A. Frasca^, M. Zusi^, F.Getman^, E. Covino^, D. Gandolfi"^; ^INAF - Osservatorio Astrofisico di 
Catania, Italy, ^INAF - Osservatorio Astronomico di Capodimonte, Napoli, Italy, ^INAF - Istituto di Astrofisica e Planetologia 
Spaziali, Roma, Italy, Landessternwarte Kdnigstuhl Heidelberg, Germany. 


We present results, recently published in [5] and based on 
FLAMES/UVES and FLAMES/GIRAFFE new spectroscopic 
observations, complemented with FORS2 and VIMOS spec¬ 
troscopy taken from the literature, of young stellar objects 
(YSOs) in the L1615/L1616 star-forming region [9]. Our Hnal 
sample consists of 54 objects with masses M* ~ 0.1 —2.3M0. 

L1615/L1616, at a distance of about 450 pc [14], is a 
small Orion outlying cometary cloud [2] in which there is 
evidence of ongoing star formation activity maybe triggered 
by the UV rediation coming from massive stars in the Orion 
OB association [13,3]. 

As a continuation of previous works performed by [3] and 
[ 8 ], mainly focused on the investigation of the star formation 
history, the relevance of the triggered scenario, and the study 
of initial mass function, here we used optical spectroscopy to 
address the accretion properties of the YSOs in the cloud. In 
particular, we considered the luminosity of Hve emission lines 
(Ha, H/3, He I A5876X, He I A6678X, and He I A7065X) 
to estimate the mass accretion rates (Macc) and investigate 
possible correlation with stellar and infrared properties. 

Our main results can he summarized as follows: 

1. The YSOs in L1615/L1616 observed with FLAMES 
show a distribution in radial velocity peaked at Vrad = 
23.2 ± 3.1 km s“^ and a mean lithium abundance of 
logn(Li) = 3.3 ± 0.3 dex. This conhrms their mem¬ 
bership to the cometary cloud and their youth. 

2. The fraction of accretors in L1615/L1516 we estimated 
through the criteria by [15] is close to 30%. This per¬ 
centage is consistent with the fraction of disks reported 
by [12] for an age of ~3 Myr. 

3. The Macc derived through secondary diagnostics (in 
the range ~ 10 “^° — 1 O“^M 0 yr“^) and the accretion 
properties of the YSOs in L1615/L1616 are similar to 
those of other star-forming regions, like Lupus or Tau¬ 
rus. This might imply that environmental conditions to 
which the cometary cloud is exposed do not influence 
the accretion evolution of the YSOs within the cloud. 

4. Different methods and evolutionary tracks used to de¬ 
rive stellar parameters and mass accretion rates intro¬ 
duce dispersion in the Macc — Mi, relation (see Fig. 1), 
as recently found by other authors. 

5. As in previous studies, strong accretors (i.e., those with 
log Macc ~ — 8.5 dex) show large excesses in the 
JHKs bands, indicative of inner optically thick disk. 
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Figure 1: Macc versus M* obtained from [4,6], [7], and [11] 
tracks, respectively, from the top to the bottom panel. Dashed 
and dotted lines represent the Macc — Mi, relation and the Icr 
deviation obtained by [1] in Lupus. The vertical arrow shows 
the position of a subluminous object we found in the sample 
and after applying the correction described in [5]. The dash- 
dotted line represents the boundary at 3 Myr due to chromo¬ 
spheric activity [10]. Mean error bars are overplotted on the 
lower-right comer of each panel. 
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MASS ACCRETION IN PRE-MAIN-SEQUENCE STARS: THE CORE STRIKES BACK 


D. Fedele^, C. F. Manara^; ^Max Planck Institut fiir Extmterrestrische Physik, Garching bei Miinchen, Germany, ^ESA/ESTEC, 
Noordwijk, The Netherlands. 


The mass accretion process in pre-main-sequence stars is 
one of the main processes that govern the dissipation timescale 
of gas in disks. Accretion produces a strong and energetic 
radiation field that regulates the gas temperature and thus it 
also controls the chemical evolution of the disk. Despite its 
importance, our understanding of the mass accretion process 
is still largely imcomplete. In particular, early measurements 
of the mass accretion rate showed a large scatter of the rate of 
mass accretion for a given stellar mass (e.g., [1]). This implies 
that viscous evolution and mass accretion are controlled not 
only by the central star. Note however that part of the scatter 
in the Macc — M* is due to observational uncertainties. 

The most direct tool of investigation of mass accretion is the 
measurement of the Balmer continuum produced at the base of 
the accretion shock. The Balmer continuum emission in excess 
to the stellar photospheric emission can be then converted to 
an accretion luminosity and finally to the mass accretion rate. 
This requires the knowledge of the stellar properties (e.g., 
mass, temperature, radius) and of the interstellar extinction. 
To reduce the observational uncertainties is important to cover 
simultaneously a wide wavelength range. This is possible 
with, e.g., X-Shooter at the VLT (ESO, Chile), that offers 
medium resolution spectroscopy {R ^ 10000 — 15000) from 
the atmospheric UV edge to the K-band in the near-infrared. 
Using X-Shooter, we performed a survey of the pre-main- 
sequence stars in the Chamaeleon I star forming region. The 
sample is made of 39 targets, including binaries with spectral 
type between early-K and late-M. Stellar properties, interstel¬ 
lar extinction and accretion luminosity are measured following 
the method described in [2]. Figure 1 shows the Macc — AT* 
diagram. For comparison, we also plot the measurements for 
the PMSs in Lupus by [3], obtained and analyzed in a similar 
fashion. As can be seen there is a systematic difference be¬ 
tween Lupus and Cha I: the Lupus PMSs show a rather tight 
correlation between mass acretion rate and stellar mass while 
the scatter is much larger for the Cha I stars. Note in particu¬ 
lar that: 1) Cha I hosts many more strongly accreting objects 
compared to Lupus and 2) the scatter in the Macc — AT* for 
Cha I seems to increase with the stellar mass. 

The comparison of the Lupus and Cha I X-Shooter dataset 
suggests that the mass accretion properties of the two regions 
are different. Observational uncertainties can be ruled out as 
both dataset are taken with the same instrumental setup and 
analyzed with the same method. What is the origin of the dif¬ 
ferent spread in the Macc — AL* diagram between Lupus and 
Cha I? The two regions have slightly different ages and because 
of viscous evolution the mass accretion rates may differ be¬ 
tween the two. We note however that the age difference is too 


small compared to the time needed for the viscous evolution to 
produce detectable variation [4]. A possible explanation can 
be the initial conditions of the cloud as suggested by [5]. They 
proposed a scenario in which the scatter in the Macc — M^, 
diagram reflects the spread of the angular momentum at the 
beginning of the collapse of the molecular cloud. Complete 
samples in both star forming regions will allow us to definitely 
confirm the existence and determine the origin of the different 
observed Macc — AT* relations. 



Figure 1: Mass accretion rates vs stellar mass determined from 
X-Shooter spectra for the sample of young stellar objects in 
the Chamaeleon I star forming region (red stars) and compared 
with the values for a sample in Lupus [3]. The best fit relation 
and its dispersion is shown for the Chamaeleon I objects with a 
red line and the red dashed lines, while the same is shown with 
black lines for the Lupus sample. Empty symbols are used for 
underluminous or non-accreting targets. Adapted from [6]. 


References. 

[1] Gullbring, E., Hartmann, L., Briceno, C., & Calvet, N. 
1998, ApJ, 492, 323 

[2] Manara, C. R, Beccari, G., Da Rio, N., et al. 2013b, A&A, 
558, A114 

[3] Alcala, J. M., Natta, A., Manara, C. R, et al. 2014, A&A, 
561, A2 

[4] Manara, C. R, PhD Thesis, LMU Munich, 2014 

[5] Dullemond, C. R, Natta, A., & Testi, L. 2006, ApJ, 645, 
L69 

[6] Manara, C. R, Redele, D., & Herczeg, G. in prep. 


18 







A STUDY OF YOUNG STELLAR OBJECTS IN THE 7 VELORUM AND CHA I CLUSTERS 


A. Frasca; INAF - Osservatorio Astrofisico di Catania, Italy. 


I present here a comparative study of the main properties 
of the young stellar objects (YSOs) belonging to 7 Velorum 
(7 Vel) and Chamaeleon I (Cha I) clusters/associations, fo¬ 
cusing on their rotation, chromospheric radiative losses, and 
accretion. All these quantities are strongly related to the stellar 
mass and the evolutionary status. 

7 Vel (age '^5-lOMyr) and Cha I (age 2Myr) are 
the first two young clusters observed in the framework of the 
Gaia-ESO Survey (GES), a large public spectroscopic survey 
aimed at the observation of 10® stars with the FLAMES 
multi-object spectrograph at the VLT/UT2 telescope [1]. 

In this work I used the fundamental parameters (effective 
temperature, surface gravity, lithium abundance, and radial ve¬ 
locity) delivered by the GES consortium in the first internal 
data release to select the members of 7 Vel and Cha I among 
the UVES and GIRAFFE spectroscopic observations. A total 
of 140 7 Vel members and 74 Cha I members were studied. 
The procedures that I developed for the GES analysis node 
of Catania Observatory to derive stellar fundamental parame¬ 
ters provided also measures of the projected rotational velocity 
(vsini) and of the veiling (r). Stellar luminosities (L*) and 
masses (M*) have been derived from spectral energy distri¬ 
butions and by the comparison with evolutionary tracks in the 
HR diagram. The equivalent widths (EWs) and the fluxes of 
the Ha and H/3 lines were measured thanks to the spectral sub¬ 
traction of low-activity and slowly rotating templates, which 
are rotationally broadened to match the v sin i of the targets. 
The Ha line was also used for identifying accreting objects, 
on the basis of its equivalent width and the width at the 10% of 
the line peak (10%tU), and for evaluating the mass accretion 
rate (Macc). 

The distribution of v sin i for the YSOs of 7 Vel is peaked 
at about 10 kms“^ and shows a tail toward faster rotators. 
This suggests the presence of both slowly-rotating stars that 
are still disk-locked or just decoupled from their disks, and 
others that have started to spin-up during the phase of contrac¬ 
tion and approach to the ZAMS. There is also some indication 
of a different v sin i distribution for the members of its two 
kinematical populations identified by [2]. A similar distribu¬ 
tion, but with a narrower peak, is observed in Cha I. The lower 
fraction of fast rotators is likely related to the younger age of 
Cha I. 

However, the most striking difference between the YSOs 
in 7 Vel and Cha I is found in the behaviour of the Ha line. 
Only a handful of stars in 7 Vel (about 4% of the members) 
display signatures of accretion, while many more accretors 
('^ 40%) were detected in the younger Cha I, where the high¬ 
est Ha fluxes are mostly due to accretion, rather than to chro¬ 


mospheric activity. Accreting and active stars occupy two 
different regions in a T^f f-Fna/Fi,oi diagram, as shown in 
Fig. 1, where a simple dividing line has been adopted to dis¬ 
tinguish the regime of Ha fluxes dominated by chromospheric 
emissiom from that where the accretion is dominant. Further 
details can be found in [3]. 

All the non-accreting stars cooler than about 4500-5000 K 
are very close to the saturation of magnetic activity, as also in¬ 
dicated by the absence of correlation between the Ha flux and 
V sin i. There is a weak tendency of a decrease of flux with in¬ 
creasing temperature for Teff> 5000 K, which is confirmed 
by the preliminary analysis of GES data of older clusters such 
as IC 4665 (age - 30 Myr) and NGC 2547 (age - 35 Myr). 

WederiveMacc values in the ranges 1O“^^-1O“®M0 yr“^ 
and 1O“^°-1O“^M0 yr“^ for 7 Vel and Cha I accretors, re¬ 
spectively. We find less scatter in the Macc — AT* relation 
derived through the Ha EWs, compared to the Ha 10%IU 
diagnostics, in agreement with other authors. However, the 
Ha 10% width represents a fast and efficient criterion to select 
accretor candidates for a deeper analysis. 

Moreover, although a wide range of Macc values is ob¬ 
served for the objects with a negligible veiling (r < 0.25), all 
the stars for which we found a veiling greater than 0.25 are 
accretors. Furthermore, the mass accretion rate displays some 
correlation with the veiling (see Fig. 2). A Spearman’s rank 
correlation analysis [4] gives a coefficient p = 0.28 with a 
significance of its deviation from zero a = 0.10, which in¬ 
dicates a marginal correlation. If we limit the analysis to the 
objects with a significant veiling (r > 0.25), the degree of 
correlation is enhanced, as indicated by the values of p = 0.62 
and a = 0.001. 

This work shows the huge potential of the GES for the 
study of rotation, activity, and accretion of YSOs in young 
clusters. This study will be extended to all the pre-main se¬ 
quence clusters that are being observed by the GES. 
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Figure 1: Ha flux (in units of the bolometric flux) versus Te// for the members of Cha I {left panel), 7 Vel {middle panel), 
and NGC 2547 {right panel) clusters. The symbol size scales with the v sin i. Filled symbols denote the accretor candidates 
{IOVcWhc > 270 kms“^) The flux values corrected for veiling by the factor (1 + r) are denoted by arrowheads. In each box, the 
dashed straight line is drawn to follow the upper envelope of the sources without accretion and is very close to the saturation of 
chromospheric activity (dotted line). 



Figure 2: Mass accretion rate from EWhc versus veiling. Diamonds and triangles represent Cha I and 7 Vel stars, respectively. 
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EXORCISM: EXOR OPTICAL-INFRARED SYSTEMATIC MONITORING 


T. Giannini; INAF - Osservatorio Astronomico di Roma, Monte Porzio Catone, Italy. 


EXORCISM is the long-term optical-infrared monitoring 
of EXors, a class of Pre-main Sequence stars undergoing inter¬ 
mittent outbursts (with timescales of tens of years and duration 
of about one year) with typical amplitudes 2-5 magnitudes in 
the visual band (see e.g. [1]). These bursts are usually inter¬ 
preted as due to episodic accretion events; nevertheless several 
fundamental aspects of the EXor phenomenon are still very 
poorly known: 1) are EXors peculiar objects or rather they 
represent a short but common and repetitive phase of the pre- 
main evolution? 2) which is the mechanism at the origin of the 
burst (disk thermal or gravitational instability, perturbation of 
the disk induced by an external body?) 3) which is the role of 
the episodic accretion in disk fragmentation and evolution? 4) 
are EXors distinct objects from FUors (pre-main sequence ob¬ 
jects with stronger bursts with longer duration) or rather they 
represent the same phenomenon in a subsequent evolutionary 
phase? 

From an observational point of view the poor knowledge 
of the EXors behavior is mainly due to a lack of long-term 
multi-wavelength monitoring programs aiming at studying the 
photometric and spectroscopic properties. Moreover, only few 
studies exist that compare photometry and/or spectroscopy of 
the outburst and quiescence phase ([2],[3],[4]) and observa¬ 
tions at high angular resolution to spatially resolve the inner 
disk. 

Our monitoring envisages to perform on monthly basis op¬ 
tical and near-infrared photometry of all the known EXors and 
candidates (around 30 objects). In addition, optical-infrared 
spectroscopy is performed roughly every year. Prompt detec¬ 
tion and follow-up of any burst is also planned. The facilities 
routinely involved in the program are : AZT24 at Campo Im- 
peratore (Italy) for near-infrared imaging and spectroscopy of 
northern objects; REM at ESO-La Silla for optical-infrared 
imaging of southern objects; LX200 at St. Petersburg Uni¬ 
versity for optical imaging. In addition, several focussed pro¬ 
grams are on-going at TNG and NOT (Canary Islands), and 
LET (Tucson) essentially to perform deep optical-infared spec¬ 
troscopy (Fig. 1). EXORCISM have so far produced a number 
of papers, to which we refer as far as the main results of our 


study are concerned ([5]-[18]). 
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Figure 1: Optical (LBT-MODS) and near-infrared (TNG-NICS) spectrum of V1118 Ori during a quiescent phase. 
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ACCRETION AND WIND PROXIES THROUGHOUT DISK EVOLUTION 


C.F. Manara; ESA/ESTEC, Noordwijk, The Netherlands. 


The evolution of protoplanetary disk and their ability to 
form planets is strictly related to the distribution and the 
amount of gas in the disk [1], This can directly be traced 
in the outer disk through observations of molecular gas emis¬ 
sion lines at various wavelengths, e.g., in the sub-mm [2], On 
the other hand, the interaction of the disk with the central star 
through accretion of matter allows us to access observationally 
to the gas content in the innermost region of the disk, at dis¬ 
tances as close to the star as ~0.1 AU [3]. Similarly, studies 
of the wind properties of young stellar objects (YSOs) help us 
in constraining the properties of the gas in the inner part of the 
disk [4] and the mechanisms that lead to the final disk disper¬ 
sal [5], Here I report on recent work carried out with various 
members of the Jets and Disks at INAF (JEDI) collaboration 
to study protoplanetary disk evolution through the analysis of 
accretion and wind properties of YSOs in various star-forming 
regions and at different stages of disk evolution. 

The data presented in this contribution have been obtained 
with the VLT/X-Shooter spectrograph, that covers simultane¬ 
ously and with medium resolution the range of wavelengths 
A <^300-2500 nm. Therefore, spectra obtained with this in¬ 
strument allow us to access the most direct proxy of accretion, 
that is the UV-excess short-ward the Balmer jump, several 
permitted emission lines (e.g., Ha, H;9, CaH, CaK, Pa/3, Br 7 
[6]) and forbidden lines (e.g., [OI]A 6300 X [4]), and also 
several photospheric absorption features (e.g., Li). We have 
recently developed an automatic procedure to determine from 
such spectra the stellar and accretion parameters of YSOs [7]. 
We find the best fit in a grid of models that includes var¬ 
ious photospheric templates, mostly from [8], a set of slab 
model to reproduce the excess emission due to accretion, and 
a reddening law to take into account the effects of interstellar 
absorption (Fig. 1). This method has already been applied to 
a large number of YSOs in the Lupus region [6], where we 
have found that the scatter of the mass accretion rates (Mace) 
at any given stellar mass (M*) is significantly smaller than in 
previous works. We have also used this method to study a 
sample of 22 transition disks (TDs) [9]. Our study has shown 
that there are (at least) some TDs with Macc as high as those of 
full disks, a finding difficult to be explained by current mod¬ 
els. Indeed, this suggest that the innermost region of their 
disk, although dust depleted, is significantly gas rich in order 
to sustain the observed Macc- Similarly, we have found that 
the winds traced by the [01] A 6300 X line scale with the accre¬ 
tion rates in the same fashion as in full disks. This line could 
trace gas from a large region in the disk [10], and we are now 
carrying out further studies with higher resolution spectra to 
understand whether this is valid also in TDs, and what kind of 
wind is this line actually tracing [11]. 

To further constrain the evolution of protoplanetary disks 
is important to study objects with different intrinsic proper¬ 
ties, such as age, M*, or the environment in which they are 
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Figure 1: Example of best-fit of X-Shooter spectra using the 
procedure described in [7]. The observed spectrum is shown 
in red, the best fit in blue, the reddened photospheric template 
in green and the slab model spectrum in cyan. 


located. With the JEDI collaboration we have studied eight 
very low-mass objects in the ~3-5 Myr old cr-Orionis cluster 
[12], 36 YSOs in the <^2-3 Myr old Lupus region [6], and 
some objects in the old (<^5-10 Myr) TWA association [13]. 
Two more samples have now been added to this collection. 
Here 1 briefly report the results on the sample of 17 very low- 
mass YSOs and brown dwarfs (BDs) located in the young 
(~1 Myr) p-Ophiucus cluster [14], while in the contribution 
by Eedele we report on the sample of 31 YSOs in the <^2-3 
Myr old Chamaeleon 1 region [15]. Studying YSOs located 
in the p-Ophiucus cluster is more challenging than studying 
objects in other older and less embedded region. Indeed, the 
extinction towards a single target in p-Ophiucus can be as high 
as Av >^10-15 mag, or even more. Therefore, none of the 
X-Shooter spectra we obtained had enough S/N in the Balmer 
continuum to perform the detailed analysis described earlier. 
We then obtained the stellar parameters (spectral type, Av) by 
finding the best match between the observed spectrum and a 
set of reddened photospheric templates, again mostly from [8] 
or from synthetic spectra. Then, we determined the accretion 
rates from the intensity of various emission lines. Eig. 2 shows 
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C.F. Manara - Accretion and wind proxies throughout disk evolution 


the values of Macc vs M* for the objects analyzed in this study 
and for the Lupus and cr-Orionis samples. The best-fit Macc- 
M* relation of the Lupus sample is also shown, and it is clear 
that most of the p-Ophiucus objects follow that same relation. 
This suggests that this relation holds in the very low-mass and 
BD range, as well. Moreover, also the values of Macc for 
p-Ophiucus targets seem to have a similar small scatter around 
the best-fit relation as the Lupus objects. It is possible that this 
scatter smaller than in previous works is related to the good 
quality of the data and the uniform and accurate methodology 
used in these works. This suggests that larger studies of un¬ 
biased samples based on similar data and analyzed with our 
methodology will really help us in determining the main prop¬ 
erties of the accretion process, and thus of the properties of 
gas in the innermost part of protoplanetary disks as a function 
of M* and age. 

We are now collecting a large number of X-Shooter spectra 
to complete the samples of YSOs with disk in the Lupus and 
Chamaeleon I regions, and we are asking for similarly com¬ 
plete data-sets in the (j-Orionis, Upper Sco, and p-Ophiucus 
regions. This will allow us to cover a large range of M* and 
age. All these spectroscopic surveys are coupled with on¬ 
going or future ALMA surveys aimed at determining the outer 
disk properties. We thus have now a detailed method available 
to determine stellar and accretion properties of YSOs, we are 
now working on determining what kind of wind is traced by 
the forbidden emission lines, and we will then be able few 
months from now to clearly constrain these processes on large 
samples of objects with well determined outer disk properties. 



log(M*/Mg) 


Figure 2: Mass accretion rates vs stellar mass for a sample of 
very low-mass stars and brown dwarfs in the p-Ophiucus clus¬ 
ter [14], together with objects in the Lupus region [6] and in 
the (j-Orionis cluster [12]. The best fit relation for the Lupus 
sample alone and its scatter is also shown. Adapted from [14]. 
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INVESTIGATING MID-INFRARED HYDROGEN LINES AS ACCRETION INDICATORS 


E. Rigliaco; Department of Physics - Institute for Astronomy, ETH Zurich, Switzerland. 


I summarize here the properties of the Hydrogen recom¬ 
bination lines observed in the mid-infrared (mid-IR) Spitzer 
spectra of disk-bearing objects [9]. In particular, I investigate 
their nature as accretion indicators, and the possibility of using 
these lines to measure accretion in embedded class-I objects 
where optical tracers cannot be used. 

A quantitative estimate of the accretion of mass from the 
circumstellar disk onto central star in young stellar objects 
(YSOs) can be measured using several indicators. In the past 
few years many works have in fact shown that many emission 
lines observed in YSOs’ spectra (mainly hydrogen recombina¬ 
tion lines) are correlated with the accretion luminosity (Lace, 
[3,7,4,5,8,2,1], among others). Even if these accretion indi¬ 
cators span over a huge range of wavelengths, from the near¬ 
ultraviolet to the near-infrared, they cannot be used to measure 
accretion in embedded objects, where the UV and optical radi¬ 
ation is absorbed by the remnant infall envelope and re-emitted 
at longer wavelengths. 

In [9] we focused on the strongest of the Humphreys Series 
lines in the mid-IR, the Hua at 12.37/im (hereafter Hl(7- 
6 )) and the Hl(9-7) line at 11.31/im, both observed with the 
InfraRed Spectrograph (IRS) on the Spitzer Space Telescope, 
aiming at unveiling their nature as accretion indicators. 

We analyzed a sample of ~100 disk-bearing objects col¬ 
lected from the Spitzer Archive in different evolutionary stages: 
from full to transitional and debris disks. We measured the Hi 
line luminosities, and the hydrogen line ratios. 

We found that full and transitional disks show a correlation 
between the Hi line luminosity and Lace [9]: 

Log{LHi(y-6 )/Lq) = 

= (0.48 ± 0.09) X Log{Lacc/LQ) - (4.68 ± 0.10) 

The same correlation cannot be tested against debris disks since 
no information about the accretion luminosity can be retrieved 
due to their more advanced evolutionary stage. Moreover, we 
found that the physical properties of the emitting gas, obtained 
through the Hl(7-6)/ Hl(9-7) line ratio are fairly well con¬ 
strained using the prediction of the theoretical models by [ 6 ] 
that replicate the conditions in the region of winds and accre¬ 
tion flows in accreting YSOs. These two results suggest that 
the Hi mid-infrared lines observed in full and transitional are 
likely accretion indicators as the other hydrogen recombination 
lines at shorter wavelengths. 

Interestingly, Hl(7-6) emission is also observed in a sam¬ 
ple of debris disks. A chromospheric origin of these lines 


seems likely, but it is not possible to definitively rule out either 
a long-lasting accretion signature nature or the presence of a 
non-resolved accreting companion. If we assume the same 
correlation in Eq. 1 can be used for debris disks, we find that 
the observed Hi lines would trace accretion luminosities as 
low as <lO~^°M 0 /yr. Optical and UV accretion diagnostics 
are not sensitive to such low mass accretion rates for Sun-like 
stars. 

The analysis of archival Spitzer spectra of disk-bearing 
objects has highlighted two main aspects that deserve to be 
better investigated in the future: 

• the presence of mid-IR hydrogen recombination lines 
in debris disks; 

• the possibility of using these line luminosities to de¬ 
termine the accretion rates in still embedded class I 
objects. 

Although in class I objects the obscuration of the protostar 
by the envelope and disk is significant (A\/=10-50 mag), it 
can still be detected at infrared wavelengths. Estimating the 
infall accretion in class I objects will certainly help to better 
understand the earlier stages of stellar evolution, and within the 
JEts and Disks INAE projects it will help to link the studies of 
pre-stellar cores and molecular jets to the accretion properties 
of more evolved class II objects. 
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EFFECTS OF ASYMMETRIC JETS ON THE DYNAMICS OF PROTOPLANETARY DISKS: STUDY OF A SIMPLE 
MODEL 


F. Bacciotti^, U. Locatelli^, R.I. Paez^, M. Volpi^; ^INAF - Osservatorio Astrofisico di Arcetri, Italy, ^Dip. di Matematica, 
Universitd di Roma “Tor Vergata”. 


I. Introduction 

Extrasolar planets can appear in very different dynamical con¬ 
figurations (see [2]). In particular, the orbits of some multi¬ 
planetary systems can show values of the eccentricities much 
larger than those in our solar system (in which the orbits are 
nearly circular). This is somewhat surprising, as large eccen¬ 
tricities could prevent the stability of such systems (because 
close encounters between planets become more likely). 

Here we investigate if the orbital properties of the planets 
can be related to the formation of these systems, and in par¬ 
ticular to bipolar collimated outflows of matter (jets), which 
are ubiquitous around Young Stellar Objects (YSOs). The jet 
physical properties can be studied using spectroscopic diag¬ 
nostic techniques (see [1]). In particular, this allows us to 
measure the flux of linear momentum (hereafter p) of the jets 
in the two lobes, that propagate perpendicularly to the disk 
midplane. The jet lobes are often asymmetric and are found to 
have a different linear momentum flux (that is Ap 7 ^ 0). This 
produces a feedback on the inner part of the disk, from where 
the jet is accelerated, because the orbital motion of the matter 
in this region is influenced by a rection force — Ap (assumed 
here to be perfectly vertical with respect to the disk midplane). 

In the present work, we want to describe the first steps of 
an ambitious project willing to link the two topics introduced 
above. In particular, we want investigate if the described jet 
feedback excites the eccentricities of the orbital motions up to 
the values observed in extrasolar systems. Said in other words, 
are the orbital properties of the currently observed exoplanets 
a legacy of the primordial formation mechanism? 

Some of these ideas have been investigated in the past by 
Namouni (see [7] and the references therein). However, in our 
model we remove some inconsistencies of that approach. In 
particular we do not consider already formed (proto)planets, 
which begin to exist close to the end of the estimated lifetime 
of jets, i.e. 10® yrs (see [3]). Moreover, differently by [7], 
we do not model the feedback — Ap as due to an asymmetric 
stellar wind, but we assume it is due to a disk wind (which has 
a much larger effect on the system). 


F = — Ap acting on the first ring as a reaction to the launch¬ 
ing mechanism producing the two lobes of the jet; (VII) all 
the forces are neglected except for F and the gravitational in¬ 
teractions between the star and each of the two rings, and the 
rings one with each other. 



Figure 1: Scheme representing the elements of the model. 

We denote with R\ and R 2 the radii of the two rings, 
while zo , zi, Z 2 , mo , m\ and m 2 are the vertical quota and 
the masses of the star and the rings, respectively. Let us em¬ 
phasize that the rings can both pulsate and move vertically, 
thus, R\, R 2 and Zj (V j = 0,1, 2) are functions of time t. 
Moreover, the evolution of each ring can be described by the 
motion of a single point belonging to it, because the position 
of all other points is obtained by suitable axial rotations start¬ 
ing from that point acting as a generic representative of the 
ring itself. This highlights the fact that our model has nine 
degrees of freedom (d.o.f.). In order to separate the motion 
of the barycenter and to suitably use our numerical integrator 
(see next section), it is convenient to adopt the Jacobi canoni¬ 
cal coordinates (vi,V2,v\, V2). Here vi locates the generic 
representative of the first ring with respect to the star, while 
V 2 is nearly equal to the astrocentric position of the generic 
representative of the second ring, while Vj the conjugate 
momentum with respect to Vj ,'tf j = 1, 2 (see [4]). There¬ 
fore, our model is fully described by the following 6 d.o.f. 
Hamiltonian: 


2. The model and the choice of its parameters 

We drastically simplify the problem of the disk dynamics, by 
making the following assumptions, which are schematically 
represented in Figure 1: (I) the (proto)star is a point mass; 
(II) both the inner part of the disk (at the jet base) and the 
external one (all the rest) are represented by two horizontal 
non-rigid circular rings of matter; (III) the whole system is 
axisymmetric with respect to the vertical axis passing through 
the star and the centers of both rings; (IV) the masses of the 
star and each ring are constant; (V) masses are homogeneously 
spread along each ring; (VI) there is a constant upward force 


H{vi,V2,Vi,V2) = A{vi,V2,Vi,V2) +B{vi,V2) , ( 1 ) 


where the main part A = ICi{vi,vi) + IC 2 {v 2 ,V 2 ) is given 
by the sum of two separated Keplerian terms, that are 




N ttj ■ Vj 
2pj-imj 


Gpj-imj 

Ikfll 


Vj = l,2, (2) 


being po = mo , pi = po + mi and p 2 — pi + m 2 ■ In 
general, the Keplerian part A is much greater than the purely 
potential term, which is 


B = -. 


Gmom2 

\\V2+^Vl\ 


Gpim2 


+ Ui -f U12 , 


( 3 ) 
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where both Vi, that is the potential energy corresponding to 
the vertical force F, and the term due to the gravitational 
interaction between the rings, i.e. 

U,, = - r 

Jo o/Rj + 2 R 1 R 2 cos-i/; + R% 

must be expressed as functions of the Jacobi variables (v\,V 2 ). 

We fix most of the parameters and the initial conditions 
as follows. According to the values commonly reported in 
the scientific literature, we set m\ + m 2 = mo/ 100 and we 
assume that ro = 0.01 AU, rt = 10 AU, = 200 AU, 
where ro and ri denote the radii of the circles bounding the 
internal (orange) jet launch region in the upper part of Figure 1, 
resp., while re is the outer radius of the disk. The initial value 
of the radius R\ (0) is set to be equal to the average distance 
of the internal part of the disk, while R2 (0) is fixed so that the 
gravitational energy due to the interaction between the first ring 
and the outer disk is exactly equal to the value assumed by U12 
in (4), when we put R\ = i?i (0) and R2 = J ?2 (0). In order to 
perform these calculations, we need to assume a disk density 
distribution law, that can be found, e.g., in [ 8 ]. Moreover, all 
the initial quotas and vertical velocities are set equal to zero, 
while the initial angular velocity of the matter in each ring is 
fixed so that its radius would constant if the other ring and 
the vertical force F are not present. Just two parameters have 
not yet been fixed: the stellar mass mo and F = —Ap. We 
compiled these values for different YSOs using from the recent 
literature (see Figure 2). We are interested in evaluating the 
effects induced by the strongest vertical force, and thus we 
focus on the case of LkHa 233, where mo = 2.9 Mq and 
||F|| 6.4 X 10“® MqAU/ yr^ (see [ 6 ]). 
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Figure 2: [|p|| (symbol: x) on a single lobe of ajet and || Ap|| 
(when data on the two lobes are available, symbol: •) as a 
function of the stellar mass. 


3. First numerical results and conclusions 

In order to numerically integrate the twelve Hamilton equations 


dvj 

dt 


Vi = 1,2, (5) 


it is convenient to use the method SBAB 3 , belonging to a 
set of symplectic integrators often used in Celestial Mechanics 
(see [5]). The tests on the preservation of the energy (as well 
as many details missing in this short note) are reported in [9]. 

Despite both rings always remain horizontal and circular, 
each generic representative of a ring moves on an orbit that 
can be instantaneously approximated with a Keplerian ellipse. 
Thus, we can numerically compute the evolution (in time) of 
the value of the eccentricity that is related to each ring. Figure 3 
shows the behavior of the eccentricity related to the first ring 
for a model, having the same values of the parameters reported 
in the previous section, but the vertical force that is increased 
by a factor 100, i.e., now ||Ft| rx 6.4 x 10“"^ M^AU/yr"^. 
In Figure 3, one can appreciate that the average eccentricity is 
nearly equal to the current value characterizing the Earth orbit. 

In conclusion, our first numerical investigations show that 
the the feedback exerted by asymmetric jets on the disks can 
effectively increase the eccentricity of orbital motions (initially 
set to be circular). However, this mechanism looks far from 
being able to explain the large eccentricities observed in some 
exoplanetary systems. Possibly this is due to a deficiency 
of our model. Our next step will be to refine the model by 
including a much larger number of rings, so to describe a 
situation closer to a disk with a continuous distribution of 
matter. 


Eccentricity of the FIRST ring 



Figure 3: Behavior of the eccentricity related to the first ring. 
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SHOCKS IN ACCRETION/EJECTION PROCESSES: OBSERVATIONS, MODELS, AND LABORATORY EXPERI¬ 
MENTS 

R. Bonito^’^; ^Dipartimento di Fisica e Chimica - Universita' di Palermo, Italy, ^INAF - Osservatorio Astronomico di Palermo, 
Italy. 


Stellar jets have been observed in a wide range of wave¬ 
lengths, from radio to optical, and X-rays, showing complex 
morphologies along the jet axis: chain of knots, bow shocks, as 
well as stationary and elongated high energy emission mainly 
detected near the driving source. We developed detailed nu¬ 
merical models to explain the knotty structure detected mainly 
in the optical [3] and the stationary X-ray emission from the 
base of stellar jets [4], both considering pulsed jets (with vari¬ 
able initial velocity of ejection) and diamond shocks formed 
in the nozzle scenario [5]. The model takes into account all 
the relevant physical processes such as the radiative losses and 
the thermal conduction which are fundamental to compare the 
model’s predictions with the observations [3]. The emission 
maps synthesized from the models can be directly compared 
with the observations performed e.g. with the Hubble Space 
Telescope in the optical and with Chandra in X-rays ([2,7]; see 
also Fig. 1). 

A new approach is also discussed which takes into account 
the effects of the magnetic field: a combination of laboratory 
experiments, numerical models, and comparison with obser¬ 


vations [1]. The opposite process of mass accretion has been 
also investigated taking advantage of the numerical magneto¬ 
hydrodynamic simulations and diagnostic tools to compare the 
spectra with the observations [6]. Future perspectives include 
the laboratory experiments of accretion shocks and a magne¬ 
tized jet with the thermal conduction included [8] from which 
both thermal and non-thermal emission can be synthesized 
from the models and compared with current multi-wavelength 
observations (e.g. radio, optical. X-rays). 
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R. Bonito - Shocks in accretion/ejection processes 



Figure 1; Emission maps synthesized from the model of a pulsed jet. Both the optical (left panels) and the X-ray emission (right 
panels) can be derived from the model to be compared with the observations. The enlargement on the right shows the stratification 
of the high energy source detected at the base of the jet and the chain of optical knots along the jet axis. 
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THE CHEMICAL INVENTORY OF HHI 


T. Giannini; INAF - Osservatorio Astronomico di Roma, Monte Porzio Catone, Italy. 


Protostellar shocks exert a profound impact in the inter¬ 
stellar medium, provoking irreversible changes in its physical 
and chemical composition and ionization degree. We have ex¬ 
ploited the X-shooter sensitivity to observe the prototype pro¬ 
tostellar shock, i.e. the Herbig-Haro object HHI (aj 2000 . 0 = 
os'* 36"* 20“27, (5 j2000 .0= -06° 45' 04''84, distance of 414 
pc). Our aims are to: 1) provide the most complete atlas of 
the chemical species involved in protostellar shocks; 2 ) obtain 
a detailed picture of the gas physical conditions that will be 
used as a benchmark for shock models; 3) understand the abil¬ 
ity of shocks to destroy dust and enhance the gaseous phase 
abundance of the refractory species. 

The target was observed during several nights between Novem¬ 
ber 2013 and January 2014, for a total integration time on- 
source of ~3.5h in the range 3 000-24 700 A (for details see 

[1]). The spectrum is extraordinary rich in emission lines, 
with more than 500 detections of fine structure ionic lines from 
species with Z < 28, HI recombination lines of the Balmer, 
Paschen, and Brackett series, several He recombination lines, 
and H 2 ro-vibrational lines with v^p < 9. The analysis of the 
HH 1 spectrum consists of the following steps; 1) derivation 
of the local extinction (Av); 2) derivation of the gas physi¬ 
cal conditions; 3) estimate of the chemical abundances of the 
atomic species. 

The local extinction was derived from ratios of lines com¬ 
ing from the same levels [1,2]. In particular, from very bright 
HI ratios of the Balmer and Paschen series, we got 0.0 < 
Av <0.8 mag. Temperature and density of the emitting gas 
were fitted for each species by applying a Non Local Ther¬ 
mal Equilibrium (NLTE) model, which considers collisional 
excitation and de-excitation and spontaneous decay as pro¬ 
cesses responsible for the line emission. The main result of 
this part of the analysis is summarized in Figure 1 ([3]). Here 
the electron temperature (black data) is plotted as a function 
of IPa„e=(IPi-i+IPi)/2, i.e. the average between the ioniza¬ 
tion potentials of the ionic states i -1 and i. This gives an idea 
of the energy range corresponding to the ionization stage i of 
the considered species. Notably, the electron temperature is 
not constant, but rather spans in the range 5 000-90 000 K. 
Moreover, we find a strong correlation between the ioniza¬ 
tion energy of the ion under consideration and the derived Te. 
Conversely, two different regimes of density can be individ¬ 
uated: a more tenuous component with rie between 1.0 10 ^ 
cm~® and 5.0 10® cm“®, probed by low-ionized atoms, and 
a denser gas component, with tie > 1.0 10 ® cm~®, probed 
by high-ionized atoms, with a tendency of this latter to a be 
also at higher temperature. To derive the fractional ionization 
Xe = Ue/nn and the fractional ionic abundance X'/X*"*"® of 
each species, we developed a ionization code that takes into 
account the following processes: direct ionization, radiative 
and dielectronic recombination and direct and inverse charge- 
exchange with hydrogen. In Figure 1 red data represent Xe vs. 


IPi. Although the data points are few, a clear stratification in 
Xe is recognizable, being the gas substantially neutral for IPi 
< 10 eV (and Te < 20 000 K) up to be completely ionized for 
IPi > 40 eV (and Te > 40 000 K). 

To derive the chemical abundances of the detected species, we 
considered ratios of bright atomic lines with the H/3 line for 
species having Te < 30000 K and the Hell 4-3 line in the 
opposite case. The observed ratios were compared with the 
intensities predicted at the derived Te, rie, and Xe. As a gen¬ 
eral result, we get abundances significantly lower both than the 
solar values ([4]) and the estimates in the Orion nebula ([5]). 
This trend is even emphasized for refractory species (C, Ca, 
Ti, Fe, Ni) which abundance is on average around 50% less 
than the solar value. This result can be interpreted in terms 
of the shock efficiency, which in HH 1 as in other protostellar 
shocks (e.g.[6],[7]), is not enough to completely destroy dust. 



Figure 1: Black/Red: Te/Xe as fitted with the excita¬ 
tion/ionization model as a function of the ionization potential. 
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CONNECTION BETWEEN JETS, WINDS, AND ACCRETION IN T TAURI STARS 


B. Nisini; INAF - Osservatorio Astronomico di Roma, Monte Porzio Catone, Italy. 


The standard paradigm for the formation of low-mass stars 
is that powerful jets and/or winds are formed during the early 
stages of stellar evolution as a mechanism to extract the ex¬ 
cess angular momentum from the system, allowing accretion 
of matter from the disk onto the star. Several mass ejection 
phenomena from active T Tauri stars (Classical T Tauri, CTT 
stars) are indeed detected using different tracers/techniques. 
The most clear evidence of mass loss are the collimated and 
extended jets observed in optical narrow band images. Such 
high velocity jets can be also easily spotted in the spectra of 
their driving source through the detection of forbidden emis¬ 
sion lines (red- or blue-shifted, with velocities >40-50 km s“^) 
of abundant species, such as [01], [SII], [Nil],[Fell], which 
span a wide range of wavelengths from UV to IR (e.g. [5]). 
This emission is usually referred to as the High Velocity Com¬ 
ponent (HVC). In addition to this HVC, spectra of T Tauri stars 
exhibit forbidden line emission at low velocity (< 40 km s“^) 
which is not spatially resolved (the so-called low velocity com¬ 
ponent, or LVC). The origin of this component is less clear; 
suggestions are that it could originate from a slow disk-wind 
or from photo-evaporating wind driven by high-energy pho¬ 
tons from the central star [6,9]. Finally, the presence of hot 
stellar winds is often testified by the prominent P-Cygni ab¬ 
sorption profile shown by many CTT in the He I l.OS^m line 
[4]. The above emissions have been separately investigated 
and discussed in several studies, however rarely the different 
mechanisms have been related together and connected with the 
accretion properties of the central star in a unified picture. 

We present here a statistical investigation on the occur¬ 
rence of different jets/winds tracers and their mutual relation¬ 
ship based on spectra acquired at VLT with the X-shooter 
instrument. These data, covering in a single exposure the 
spectral range from UV to near-IR, allow us to simultaneously 
observe the main tracers of both high velocity jets, winds and 
accretion signatures. The so far considered sample is com¬ 
posed by CTTs in three star forming regions, namely Lupus, 
Chameleon and cr Ori. The X-shooter spectra of Lupus and 
(T Ori (44 sources) are presented in [1] and [8], where their 
accretion and stellar properties have been derived. The prop¬ 
erties of their LVC emission has been instead discussed in [7]. 
The Chameleon sample comprises 16 sources selected from 
the objects investigated in [2]. All together, the global sample 
of 60 sources spans a mass range between 0.03 and 2 M©, and 
mass accretion rates from 2.5 10“^^ to 1 10“^ MQ/yr. 

The most common forbidden line in the spectra is the 
[OI]6300A detected in 87% of the sources. This emission 
is more often related to a LVC (peak velocity < 40 kms“^), 
while only 32% of the sources present an HVC (peak velocity 
> 40 km s“^). Thus only ~ 1/3 of the considered sources do 
show a high velocity jet, although this number could be higher 
due to orientation effects. We have investigated if the sub¬ 
sample of sources with HVC emission present peculiarities 
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Figure 1: Correlation between [OI]6300A luminosity and 
mass accretion rate for sources of our sample in the Lupus 
(black), a Ori (green) and Chameleon (red) star forming re¬ 
gions. The line luminosity has been separately plotted for the 
LVC (lower panel) and HVC (upper panel). The correspond¬ 
ing correlation coefficient is given in each panel. For sources 
where the [OI] line is not detected in one of the two compo¬ 
nents, its luminosity has been artificially settled to 10“® L©. 


with respect to the rest of the sources and if their LVC is 
influenced in some way by the jet presence. To do that, we 
have performed a series of correlations between the extinction 
corrected [OI] line luminosity, separately computed for the 
HVC and LVC through gaussian line fitting, and accretion 
and stellar parameters. As an example, in Fig. 1 we show 
the correlation between the line luminosity and the source 
mass accretion rate (derived in [1,9,2]). The line luminosity 
correlates with M for both the LVC and HVC, however there 
is a significant number of high accretors that do not have any 
detected HVC. In practice, only about half of sources with M 
> 10“® M./yr show evidence of a high velocity jet. Similar 
correlations have been found between [OI] luminosity and 
other stellar parameters, such as mass and luminosity, without 
finding any segregation of sources without HVC in specific 
ranges of parameter space. Similarly, there is no difference in 
line luminosity or kinematics (e.g. velocity peak or line width) 
of the LVC in sources with and without the HVC. 

In the assumption that all protostars develop jets during 
the earliest phases of their formation, the above findings sug¬ 
gest that jets and disk accretion do not evolve with the same 
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Figure 2: GIANO continuum-normalized observations of the 
Hel 1.08/rm and Pa 7 lines of XZ Tau taken in sept. 2014 
(black line) compared with spectra of the same lines taken dur¬ 
ing 2002 at Keck II [4]. The recent increase in the Pa 7 EW, 
likely caused by an accretion event, is followed by an increase 
in the Hel emission component as well as by a widening of the 
blue-shifted absorption component, suggesting that the stellar 
wind is directly driven by accretion. 

timescales, but high velocity jets are dissipated while the star 
is still accreting. In addition, high velocity jets and slow winds 
seem to be unrelated, since the properties of slow winds are 
not influenced by the jet presence. 

We have also searched for the presence of P-Cygni ab¬ 


sorption in the Hel 1.08/im line of our sample, in order to 
correlate the hot stellar winds traced by this line with the slow 
and high velocity winds/jets probed by the [01] LVC and HVC. 
Only 19/60 sources have a detectable P-Cygni absorption, a 
result that can be biased by the limited spectral resolution (R~ 
8000) that prevents the detection of the weak features. For 
these sources, we find a weak correlation between [OI] lumi¬ 
nosity and the blue-shifted absorption EW, for both the HVC 
and LVC. Similarly, such EW only weakly correlates also with 
the accretion luminosity of the sources (correlation coefficient 
0.5). Such correlations would need to be better confirmed on 
a larger sample than analysed here. 

A direct evidence that the depth of the Hel 1.08/rm blue- 
shifted absorption is correlated with accretion comes from 
high resolution near-IR spectra acquired with the GIANO in¬ 
strument (R ~ 50 000) at TNG on the EXor variable source XZ 
Tau. Fig. 2 shows the recently obtained GIANO spectra in the 
Hel and Pa 7 lines, as compared with the spectra of the same 
lines acquired in a different epoch by [4]. The recent spectra 
show a large increase in the Pa 7 equivalent width accompanied 
by an increase of the Hel emission component. At the same 
time the depth and and maximum velocity of the blue-shifted 
absorption component are much larger than before. These 
observations suggest that the source had an accretion burst 
(testified by the Pa 7 increase) accompanied by the ejection of 
a stronger stellar wind (deeper absorption in a larger velocity 
range). 
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DISKS AND JETS WITH ALMA: SCULPTING THE BIRTHPLACE OF EXOPLANETS 


Linda Podio; INAF - Osservatorio Astrofisico di Arcetri, Firenze, Italy. 


The first step to comprehend the formation of our solar sys¬ 
tem as well as of extra-solar planetary systems is to study the 
evolution of the planet formation site: the circumstellar disk. 
To understand how disks lose their mass and form protoplanets 
is a front line in astrophysics. These processes are believed to 
be regulated by the simultaneous effects of the mass accretion 
onto the star and the ejection of matter from the stellar-disk 
system. However, the interface between jets and disks has 
been poorly investigated to date because of both the very small 
angular scales involved and the occurrence of numerous other 
kinematic components in the circumstellar region (cavities of 
swept-up material, infalling envelope, static ambient cloud). 
The Atacama Large Millimeter Array (ALMA) with its un¬ 
precedented combination of angular resolution and sensitivity 
offers us for the first time the opportunity to simultaneously 
probe the disk and the jet in embedded protostars. 

The HH 212 protostellar system in Orion (at 450 pc) is an 
ideal laboratory to investigate the interplay of infall, outflow 
and rotation in the earliest evolutionary phases of the star¬ 
forming process. HH 212 is low-mass Class 0 source driving 
a symmetric and bipolar jet extensively observed in typical 
molecular tracers such as H 2 , SiO, and CO (e.g., [14]). Obser¬ 
vations of the HH 212 protostar in SO Qg — 87 and lOn — lOio 
lines have been performed using ALMA-band 7 extended con¬ 
figuration (Cycle 0). The obtained data-cubes have a spectral 
resolution of 488 kHz (0.42 — 0.43 km s“^), a typical beam 
FWHM of 0.65" X 0.47" at PA~ 49°, and an rms noise of 
~ 3 — 4 mjy/beam in the 0.43 km s“^ channel. 

The channel maps presented in Figure 1 show that SO 
Og — 87 traces different components in the inner ~ 1300 AU 
of the HH 212 protostellar system: (i) the X-shaped cavity 
walls at systemic velocity; (ii) the poorly collimated outflow 
(width ~ 135 AU) swept out by the jet at low velocities; (iii) 
the collimated {Rjet ~ 90 AU) and fast (Vjet ~ 100 — 200 
km s“^) molecular jet at intermediate and high velocities. 
Conversely, the emission in the high-excitation SO lOn — lOio 
line presented in Figure 2 shows a compact structure with 
the blue- and the red-shifted peaks displaced roughly on the 
equatorial plane, which is consistent with an origin in the 90 
AU disk rotating around a ~ 0.2 — 0.3 Mq star previously 
detected by [7] and [2]. 

The ALMA observations show that SO emission at high- 
velocity can be used as an effective probe of the molecular jet, 
similarly to a standard tracer such as SiO (see also [8,13,3]). 
The high velocity molecular jet is dense (~ 10® — 10® cm~®) 
and transports a mass loss rate > 0.2 — 2 x 10 “® Mq yr“^, 
indicating a high ejection efficiency (> 0.03 — 0.3, for an 


infall rate of ~ 6 x 10“® M© yr“^, [9]). The SO abundance 
in the jet (~ 10 “^ — 10 “®) is similar to what estimated in 
other Class 0 outflows and jets (e.g., [1,13,8]) and indicates 
that 1% to 40% of the elemental S has been converted into SO 
in the shocks along the jet (e.g., [ 6 ]) and/or at the base of the 
wind due to ambipolar diffusion [ 10 ]. 

The disk mass derived from C^^O 3-2 is > 0.002 — 0.013 
M 0 , which is consistent with what inferred from HCO^ by [7]. 
The SO abundance (~ 10“® — 10“^) is 3 — 4 orders of mag¬ 
nitude larger than what estimated in evolved protoplanetary 
disks (e.g., [5]). Such an SO enhancement may be produced 
in the accretion shock occuring at the interface between the 
envelope and the disk [7,12], or in spiral shocks if the disk 
is partly gravitationally unstable (e.g., [4]). Higher angular 
resolution and sensitivity observations are required to address 
the exact origin of the high-excitation SO emission. 
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Figure 1: Channel maps of SO 9g — 87 at systemic, low (LV), intermediate (IV), and high (HV) velocities (panels from left to 
right). The green, blue, and red contours trace the emission at systemic velocity and over the labeled blue- and red-shifted velocity 
intervals {\V — I4ys|). The lowest contour of the emission at systemic velocity is shown in all panels. The first contour is at 5 (j for 
SO and 3 (t for SO 2 with steps of Sa (Vsys) and 20cr (LV, IV, HV) for SO and 3a for SO 2 (Icr ~ 3-4 mJy/beam/0.43 km s“^). The 
tilted black cross is centred at the source position and indicates the jet and the disk PA (PAjet = 22°, PAdisk = 112°). The ellipse 
in the bottom-left comer shows the beam HPBW of the line emission maps (0.65” x 0.47”, PA = 49°). The triangles indicate the 
emission peaks along the blue and red lobes. Adapted from [11]. 



Figure 2: Blue- and red- shifted emission of SO lOn — lOio at low velocities (\V — Vsys\ < 1.9 km s“^, left panel) and in the 
higher velocity range where the emission is consistent with a rotating disk (HV^ia*, : 1.9 < \V — Vsysl < 3.5 km s“^, middle 
panel) is compared with C^^O 3 — 2 emission in the range (right panel, [2]). The tilted black cross indicates the jet (PA = 

22°) and the disk (PA = 112°) direction. The ellipse in the bottom-left corner shows the beam HPBW of the line emission maps 
(0.65” X 0 . 47 ”, PA = 49°). First contour at 5 a with steps of 5 a for C^^O 3 — 2 and at 3 a with steps of Icr for SO lOn — lOio. 
Adapted from [11]. 
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I. Introduction 

Current planet formation models propose two possible mecha¬ 
nisms to explain the variety and the large number of planetary 
systems observed so far (for a review of exoplanets proper¬ 
ties, see [1]): the so-called core accretion scenario, in which 
planets are formed through a process of slow accumulation 
of solids into planetesimals and rocky cores followed by gas 
accretion; or the gravitational instability scenario, in which 
planetary mass objects are formed by the exponential growth 
of instabilities in the disk. Furthermore, the Kepler mission 
and recent meteoritic studies suggest that the core accretion 
scenario is playing an important role for many systems, sup¬ 
porting the view that small rocky bodies are assembled by 
collecting grains and pebbles that grow throughout the life of 
the protoplanetary disk [2], 

2. Assessing the dust grain size 

The first step of planet formation in the core accretion scenario 
is the growth of dust to form large aggregates and eventually 
planetesimals (see [3] for a recent review). This first phase 
of growth from ISM-sized particles to mm-sized pebbles is 
directly observable in the disk midplane by measuring the 
dust opacity as a function of wavelength in the millimetre 
domain (e.g. [4]). 

At (sub-)mm wavelength the disk emission is mostly op¬ 
tically thin and can be written as 

F,, oc S,/(r)(l — e“'^) , (1) 

where S and T are respectively the dust opacity, the dust 
surface density and the local disk temperature. Assuming a 
dust grain size distribution n(a) oc a“'*foramin < a < Umax, 
it is possible to show that the dust spectral index /3 (defined as 
Ki, oc depends on Umax only (provided that Umiii < 1/rm, 
see [5]). According to Eq. (1), the slope amm of the (sub- 
)mm SED (F^ oc is related to the dust opacity spectral 

index /3 simply as Omm ~ j3 + 2. As a conclusion, from 
(sub-)mm multi-wavelength observations is possible to directly 
estimate the level of grain growth (umax) in the disk midplane 
through the measure of Omm. Typical values of ai_ 3 mm 
3.7 are found for the relatively small ISM dust grains (Omax ~ 
0.25/im, see [5]), while [6] showed that dust grains with sizes 
Umax > 1mm lead to /3 < 1 (ai_ 3 mm < 3), regardless of the 
dust composition, porosity and size distribution. 

3. A multi-wavelength fitting tool 

In order to constrain the dust grain size distribution in a pro¬ 
toplanetary disk, we have developed a fitting tool that aims 
at fitting multi-wavelength (sub-)mm interferometric observa¬ 
tions with a self-consistent modeling of the dust opacity and 


the disk structure. The tool fits directly the observed visibili¬ 
ties in order to avoid possible nonlinear artifacts introduced by 
the CLEAN algorithm. Fitting multi-wavelength observations 
allows to break the degeneracy in the terms on the right hand 
side of Eq. (1), namely k^(R), S(J?) and T{R). 

The disk emission is modeled with a classical two lay¬ 
ers disk model (see [7]) assuming a parametrization for the disk 
surface density E(i?) = Eq (R/Rq) ’ exp [— (R/Rc)^ 'j, 
where Eq is a constant, Ro = 40 AU and Rc is the spatial 
scale of the exponential cutoff. We also parametrize the 
maximum dust grain size in the midplane as amax(F) = 
uomax (wheteas Umin is kept fixed throughout 

the disk). Therefore the disk model has five free parameters: 

y, Eq, Rc, UOmax, bmax. 

To perform the fit we adopt a bayesian approach, using 
an affine invariant Markov Chain Monte Carlo (MCMC) al¬ 
gorithm, highly efficient for multi-dimensional problems and 
designed for highly parallel computations (we use the imple¬ 
mentation of the emcee Python package, see [8]). The resulting 
MCMC provides estimates of the posterior distribution func¬ 
tion in the five-dimensional space of parameters, out of which 
we can compute the marginalized distributions (to infer the 
values for each parameter) and the bivariate distributions (to 
investigate correlations between the parameters). 

We have successfully applied this fitting tool on ALMA, 
CARMA, SMA and EVLA observations of disks from the 
Disks@EVLA sample (PI: C.Chandler, NRAO), mainly from 
the Taurus star-forming region. For each disk, we have fit¬ 
ted simultaneously four to six observations at wavelengths 
ranging from 0.8 mm to 1 cm. In almost all the cases we 
have reported j3{R) profiles that increase with radius, from 
/3 <C 1 at J? < 40 AU to /3 > 1.5 in the outer disk, and ac¬ 
cordingly Umax {R) decreasing with radius, with Umax ~ 1 cm 
at J? < 40AUand6max ~ — 1. These results confirm theoret¬ 
ical expectations from dust evolutionary models that predict 
the smaller dust particles to be diffused throughout the disk 
and the larger pebbles to be driven to the inner disk by the gas 
drag. This analysis is still statistically limited by the number of 
disks that have been analysed (less than ten), however a future 
analysis on larger samples of disks will allow to investigate 
possible correlations between the dust processing level and the 
stellar parameters. 
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